Nuclear Lattice Simulations with EFT 



Dean Lee a 

a Department of Physics, North Carolina State University, Raleigh, NC, 27603 

This proceedings article is a summary of results from work done in collaboration with Bugra Borasoy and 
Thomas Schaefer. We study nuclear and neutron matter by combining chiral effective field theory with non- 
perturbative lattice methods. We present results for hot neutron matter at temperatures 20 to 40 MeV and 
densities below twice nuclear matter density. 



1. Introduction 

Weinberg extended effective field theory meth- 
ods to the nucleon-nucleon interaction 1;, and 
over the last several years effective field theory 
methods have been applied successfully to the two 
and three-nucleon system. Our aim and the goal 
of the Nuclear Lattice Collaboration as a whole 
2 3 is to extend effective field theory methods 
to the nuclear many-body problem. For this pur- 
pose we investigate the many-body physics of low- 
energy nucleons and pions on the lattice. Our 
starting point is the same as that of Weinberg. We 
begin with the most general local Lagrangian in- 
volving pions and low-energy nucleons consistent 
with translational invariance, isospin symmetry, 
and spontaneously broken chiral symmetry. This 
yields an infinite set of possible interaction terms 
with increasing numbers of derivatives and/or nu- 
cleon fields. Degrees of freedom associated with 
anti-nucleons, heavier mesons such as the p, and 
heavier baryons such as the A, are integrated out. 
The contribution of these particles appear as co- 
efficients of local terms in our pion-nucleon La- 
grangian. We also integrate out nucleons with 
momenta greater than ira^ 1 , where a is the lat- 
tice spacing. 

The operator coefficients in our effec- 
tive Lagrangian are determined by fitting to 
experimentally-measured few-body nucleon scat- 
tering data at zero temperature. The depen- 
dence on the lattice spacing is described by 
the renormalization group and can absorbed by 
renormalizing operator coefficients. In this way 



we construct a realistic simulation of many-body 
nuclear phenomena with no free parameters. In 
our discussion we present results for hot neutron 
matter at temperatures 20 to 40 MeV and den- 
sities below twice nuclear matter density. This 
proceedings article is a summary of results from 
work done in collaboration with Bugra Borasoy 
and Thomas Schaefer 0]. 

2. Non-perturbative effective field theory 

Effective field theory provides a systematic 
method to compute physical observables order 
by order in the small parameter Q/A x , where 
A x is the chiral symmetry breaking scale and 
Q = (q^^, • • .)• Here, q is a small external mo- 
mentum and m n is the mass of the pion. The 
simplest processes are those that involve only pi- 
ons and external fields. In this case the effective 
field theory is perturbative. At any order in Q 
there are only a finite number of diagrams that 
have to be included. At lowest order these are tree 
diagrams with the leading order interaction. At 
higher order, diagrams with more loops or higher 
order terms in the interaction have to be taken 
into account. 

Weinberg showed HO that the simple dia- 
grammatic expansion for nucleon-nucleon scat- 
tering is spoiled by infrared divergences. He 
suggested performing an expansion of the two- 
particle irreducible kernel and then iterating the 
kernel to all orders to produce the scattering 
Green's function. It was later pointed out that 
a possible difficulty arises because at any order 



2 



in Q an infinite number of diagrams is summed, 
and it is not clear that all the cutoff dependence 
at that order can be absorbed into counterterms 
that are present at that order [5]. This problem 
does indeed arise if one considers nucleon-nucleon 
scattering in the 1 Sq channel |7], but in practice 
the cutoff dependence appears to be very weak 

In this work we go one step further and consider 
the nuclear many-body problem. We expand the 
terms in our action order by order, S — Sg + Si + 
S2 + • • • • At order k in the chiral expansion, we 
calculate observables by evaluating the functional 
integral 
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We will refer to this approach as non- 
perturbative effective field theory. The inter- 
actions at chiral order k or less are iterated to 
arbitrary loop order. The functional integral is 
computed non-perturbatively by putting the pion 
and nucleon fields on the lattice and using Monte 
Carlo sampling. Since the number of diagrams 
at a given chiral order grows exponentially with 
the number of nucleons, a non-perturbative tech- 
nique such as this is needed for systems with more 
than just a few nucleons. 

Computing the path integral corresponds to 
summing an infinite set of diagrams. As in the 
case of iterating the two-particle irreducible ker- 
nel to determine the full two-nucleon Green's 
function, it is not clear that the cutoff depen- 
dence at a given order in the low energy expan- 
sion can be absorbed into a finite number of co- 
efficients in the action. In practice we will there- 
fore restrict ourselves to lattice cutoffs that satisfy 
ira~ l < A x . 

3. Lowest order interactions 

We let N represent the nucleon fields. We use 
Tt to represent Pauli matrices acting in isospin 
space, and we use a to represent Pauli matrices 
acting in spin space. Pion fields are notated as 7Tj. 
We denote the pion decay constant as FP hys ps 
183 MeV and let D = 1 + 7r?/F^. The lowest 



order Lagrange density for low-energy pions and 
nucleons is given by [H], 
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gA is the nucleon axial coupling, and ey& is 
the Levi-Civita symbol. The chemical potential 
\x controls the nucleon density and [i will be set 
very close to tun. At next order we have terms 
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We will include this kinetic energy term from C^- 1 ' 
in our lowest-order Lagrange density so that we 
get the usual free nucleon propagator. In this 
study we limit ourselves to the interactions of 
neutrons and neutral pions. 

4. Lattice Schrodinger equation and phase 
shifts 

We adjust C, the coefficient of the NNNN 
contact interaction, so that the NN s-wave scat- 
tering length matches the experimental value (see 
for example JO]). In order to calculate the phase 
shifts, we will solve the lattice Schrodinger equa- 
tion for the two-neutron system and observe the 
asymptotic form of the scattering wavefunctions. 

With the lowest order two-nucleon potential 
we can construct a matrix representation for 
the Hamiltonian in the two-neutron sector and 
solve the time-independent lattice Schrodinger 
equation. At this stage one could implement 
Liischer's formula for measuring phase shifts in a 
cubical periodic box |12j . However in our case 
we can construct the eigenvectors explicitly using 
Lanczos iteration, and so we find it more straight- 
forward and accurate to read the phase shifts di- 
rectly from the asymptotic forms of the s-wave 
scattering states. We use the lattice Schrodinger 
technique to tune the coupling C to reproduce the 
large scattering length that is observed in nature. 
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Figure 1. Energy per neutron in MeV for temper- 
atures 25 MeV and 37.5 MeV and different lattice 
volumes. The inverse lattice spacing is a -1 = 150 
MeV and a t = 1.0. 



5. Results 

We have generated simulation results for a -1 = 
150 MeV; a t = ^ = 1.0; temperatures = 
25.0 MeV and 37?5 MeV; and lattice sizes 3 3 , 4 3 , 
and 5 3 . Half- filling at this lattice spacing oc- 
curs at p = 2.64p p ^ vs , where p p ^ ys is the normal 
nuclear density of about 0.17 nucleons per fm 3 . 
The calculations were performed using the zone 
determinant method |13) using the second order 
approximation with zones consisting of a single 
spatial point. By calculating the exact determi- 
nants of some generated matrix configurations we 
estimate the systemic error for the zone expansion 
to be about < 0.1% for T phys = 37.5 MeV and 
< 0.5% for T phys = 25.0 MeV. 

We have dealt with the complex action by com- 
puting the phase as an observable, and for the 
various simulations presented here we found an 
average phase of about ~ 0.95 and so this did not 
present a significant computational problem. 

In Fig. ^we show the energy per neutron as a 
function of neutron density. Our results indicate 
a rather flat function for the energy per neutron 
as a function of density. 

The flatness of our energy per neutron curves at 
these temperatures are intriguing and hopefully 



will be checked by others in the near future. The 
results of |3J also see a flattening of the energy 
per neutron curve with increasing temperatures. 
We can also compare with variational calculations 
|14| and recent quantum Monte Carlo results from 
|15j . They observe a significant flattening of the 
energy per neutron curve due to iteractions even 
at zero temperature. 
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